Functions of linear operators: Parameter differentiation 
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g ■ Abstract 

P^ . We derive a useful expression for the matrix elements [ "''%^''^^ ]i j of the deriva- 

, ^, . tive of a function /[j4(t)] of a diagonalizable linear operator A{t) with respect 



\df[A{t)]^ 



to the parameter f. The function /[^(i)] is supposed to be an operator acting 
on the same space as the operator A{t). We use the basis which diagonal- 



dt 



Ml /(^i) - /(^O In 



ij ~ ^dtiij \j -X 



\^ . izes A{t), i.e., Ai j = Aj 5ij, and obtain 

o, 

0^ , addition to this, we show that further elaboration on the (not necessarily 



simple) integral expressions given by Wilcox 1967 (who basically considered 
/[A(t)] of the exponential type) and generalized by Rajagopal 1998 (who ex- 

Ch , tended Wilcox results by considering /[^(t)] of the g-exponential type where 

expg(x) = [1 + (1 — q)x]^'^^^'^' with q £ TZ; hence, exp]^(x) = exp(a;)) yields 
^ , this same expression. Some of the lemmas first established by the above an- 

^ thors are easily recovered. 
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The quite ubiquitous necessity to have an expression of the parameter derivative of 
a general function of a hnear operator has stimulated many authors to develop various 
algorithms focusing on this problem. Special effort has addressed the calculation of the 
derivative of exponential functions of such operator in terms of the derivative of the exponent. 
Most of the results that have been achieved provide the derivative as an integral expressionlllfl, 
or by means of an integral representation^ of the function in the complex plane. Such elegant 
procedure is very useful to prove mathematical properties but, in many cases, to calculate 
the remaining integrals is not an easy task. Therefore, it is convenient to have explicit 
expressions of the matrices elements of the relevant operators. The choice of an appropriate 
basis set considerably simplifies the calculations; indeed, once these elements are obtained, 
the expression in a different basis can readily be calculated. This is the purpose of the 
present work. 

Let A{t) and A'(t) = dA(t)/dt respectively be a general diagonalizable linear operator 
and its t-parameter derivative [t does not necessarily denote time); in general they do no 
commute. Also, let /[A(t)] be an arbitrary, well behaved, function of the operator. We want 
to calculate L^ '' . For convenience, we shall express the matrices of all operators in the 
basis which diagonalizes A(t). More specifically, we have that 

[Ait)],,=X,S,, (1) 

where the {Aj} are the eigenvalues of A{t). For instance, in the usual case of perturbation 
theory in Quantum Mechanics we have 

H{t) = Ho + tHi (2) 

with Hq = H(t = 0) the unperturbed diagonal Hamiltonian, and Hi the perturbation 
of the system. Obviously for this case we have 

We shall start by considering a simple case, namely f[A] = A^ with r G 7^. The derivative 
of the power r of the linear operator A{t) is by definition 



dt s^o s ^ ' 

In order to find an expression for [A(t + s)Y it is convenient to find tlie basis in which 
A{t + s) is diagonal. More precisely, we have 

A{t + s)C = C A (5) 

where A is the diagonal matrix formed with the eigenvalues oi A{t + s) and C is the matrix 
whose columns are the eigenvectors also of the matrix A{t + s), hence 

A{t + s) = A{t) + A'{t) s + ... = A + A'{t) s + ... 

C = / + CWs + (6) 

A = A + A(')s + .. 

where we have kept terms up to order s, I is the matrix unity and A denotes the diagonalized 
matrix A{t). Identification of the coefficients of s implies 

A' + \ C(i) = C(i) A + A(i) (7) 

Then, by taking matrix elements in the last equation, we obtain 

A'^^ + KCl^ = Cl^X, [1^3) (8) 

hence 

Once we know the expression for C as a function of s we can obtain an expression for 

\A(t + s)]'', namely 

\A{t + s)X = C S!' C-^ (9) 

Then, by keeping up to terms linear in s, we can write 



A'{t + s)r^[I + s C(^)] [A + s X^^^Y [I + s C^^Y^ 



(10) 



and, since we have that [I + s C^^^] ^ = I — s C^^^ +9{s'^), we can replace this in Eq.@ and 
obtain 



d[A{t)Y _ [I + s C(^] [A + s X^^^Y [I - s g(i)] - A^' 
dt s~,o s 



C(i)A^'-A^'C«+rA(^-i)A« 



(11) 



The matrix elements obtained from the last equation are 



rO[/i{l)\ 1 _ w^_i^ (1) _ w^_i^ , 




[ g^ \ii r^^\^> \ '^I'Vj A^, 






-A, 


lave used that C\ ] = '{ . 

1' J Aj—Ai 





(^ ^ J] 



(12) 



The straightforward application of the above procedure to an arbitrary well behaved 
function of the linear operator /[y4(t)] leads us to the following expressions: 



dim)] 



dt 



^i " dX 



X=\i 



(13) 



dflAit)] 



dt 



* J ■' Xi — A,- 



An alternative proof of these same results can be obtained as followsa. Considering 
that any function /[yl(t)] commutes with the operator A{t), i.e. [ A{t), f{A) ] = 0, and 
then taking the derivative with respect to the parameter t of this commutator, and finally 
exhibiting the matrix elements, we also arrive to Eqs. (13). These two equations constitute 
the main result of this work. In what follows, and in order to exhibit their usefulness, we 
will use these expressions to prove some of the theorems given in Wilcox's and Rajagopal's 
papers. 

As a first application let us consider the formula (2.1) of Wilcoxlil, namely 



dexp[A(t)] /-i 

^^ ^ '^ = / dx exp[A{l - x)] 

at Jo 



dA{t) 
dt 



exp(a;A). 



(14) 



The matrix elements in the basis that diagonahzes the operator A{t) are given by 



(9exp[A(t)] 



dt 



« J 



dx exp[(l — x)\i 



dA{t) 



dt 



« J 



exp(a;Aj) 



dA{t) 



dt 



« J 



exp(Aj) / dx exp[x(Aj — Aj 



dA{t)i exp(Aj) — exp(Ai 
dt J « i 



Aj — Ai 



(15) 



, exp(Aj) - exp(Ai 



Aj — Aj 



which is the same expression obtained from our general formula, Eq. (13). 

As a second application of our results, let us focus on the proof of the following lemma: 



ii[B,A{t)] = f then [BJ{A)] 



dfjA) 

at ■ 



Proof: from the first equation we have 



^ij^j ~^ AiBij — Bij[Aj — Ai) — A^ 



«j 



(16) 



Taking the matrix elements of the first member of the second equation we get 



B 



« J 



/(A,) - /(A.) 



A' 



« J 



(A, - A.) 



/(A,) - /(A.) 



df{A) 



dt 



*j 



(17) 



that is what we wanted to prove. 

As a third application, let us now assume that the operator A depends on two parameters, 
namely t and u. We will prove the following theorem 



Tr 



h{A) 



df{A) dg{A) 



Tr 



h{A) 



dg{A) df{A) 



dt du 



dt du 

where /i, / and g are arbitrary functions of the operator A{t, u). 
Proof: the first term of the last equation is 



(18) 



EE MA. 



« ] 



UllK^^i) 



« 3 






dt 



du 



3 « 



A', 



* J 



(A, - A.) 



/(A,) - /(A.) 



A, 



J « 



(A. - A,) 



^(Ai) - g{\j) 



EE MA. 



« J 



A' 



« J 



(A, - A.) 



^(Aj) - ^(Ai) 



A 



J « 



EEMA 






« J 



9t 



« J 



9m 



J « 



(19) 



2^^ 



/.(A) 



dg{A) df{A) 



dt du 



''-f andi^fO 



as we wanted to prove. (We have used the notation A' 

Let us finally focus on our fourth and last application. As mentioned before, RajagopalBI 
has generalized Wilcox's expression in which the exponential function of the operator is 
replaced by a monomial fractional power of the form 



QT{t,(3) = 1 - {l-q)(3A{t) 



l-q 



expq[-(3A{t)\ 



(20) 



Such q-exponential expressions arise naturally when we consider ensembles in the context of 
the recently introduced nonextensive thermostatisticsH (for a recent review on the subject 



se 



Before we show that further elaboration on Rajagopal's expression^ gives the same result 
as the one obtained from our Eq. (13), we need the following identity 



d_ 
dx 



1 — ax 



1 — bx 



c {h — a 



[1 — ax] 



c-l 



[1 - hx] =+1 



(21) 



Defining again A' = J. , Rajagopal's equality, besides a factor P that is lacking in his 



at ' 
equation and also in Wilcox's paper, can be written in the form 

dQT{t,P) 
dt 

= -q(3 f dxQT{t,(3)[QT{t,x)]-^[l - (1 - q)xA{t)r^A'\l - (1 - g)xA(t)l"'Qr(t, x) (22) 

JO J L J 



The matrix elements of the last equality give us 



dQT{t,(3) 



dt 



qP dx 1-{1- q)p\i '"' 1 - (1 - q)x\. 



Using the identity given in Eq. (|2T|) we obtain 



1-9 



-1 



A^ h - {1 - q)xX, 



13 



1-9 



-1 



r9QT(t,/?) 



dt 



« J 



(23) 



-q(3l-{l- q)(3\ 



'-'A'^j / dX' 



1 — (1 — q)x\j 



-r — 1 

1-9 



1 - (1 -q)xXi 



T^+l 



A 



« J 



Xj — Xi 



[l-(l-q)pX,]—-[l-il-q)pX.] — 



(24) 



that is the expression that we would obtain from Eq.(13) applied to the function, QT{t,(3), 
defined in Eq.(20). 

Summarizing, we have given a general expression, Eq.(13), for the parameter differen- 
tiation of a generic function f{A) of a diagonalizable linear operator A(t) that depends on 
a parameter t. This expression enables the consistent recovering, as particular instances, 
of various useful Wilcox formulaeu, as well as their recent generalization by Rajagopala. 
We believe, consequently, that the present expression can be useful for quantum calcula- 
tions in standard or deformed Quantum Mechanics, extensive or nonextensive Statistical 
Mechanics^, among other areas of physical interest. 

We are deeply indebted to A.K. Rajagopal for a critical reading of our manuscript. 
This resulted in the beautiful alternative proof presented just after Eqs. (13), which we 
are including here thanks to his generous authorization. Partial support from CNPq and 
PRONEX/FINEP (Brazilian agencies) are acknowledged as well. One of us (D. P.) also 
acknowledges warm hospitality at the Centro Brasileiro de Pesquisas Fisicas. 
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